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The PHENIX experiment at RHIC has measured the invariant cross section for u-meson produc- 
tion at mid-rapidity in the transverse momentum range 2.5 < pt < 9.25 GeV/c in p + p and d+Au 
collisions at y / sjvjv = 200 GeV. Measurements in two decay channels (uj — > tt tt + it~ and uj — * 7T°7) 
yield consistent results, and the reconstructed uj mass agrees with the accepted value within the pt 
range of the measurements. The uj/n {> ratio is found to be 0.85 ± 0.05 stat ± 0.09 sys in p + p and 
0.94±0.08 stat ±0.12 sys in d+Au collisions, independent of p T - The nuclear modification factor R% A 
is 1.03 ± 0.12 stat ± 0.21 sys and 0.83 ± 0.21 stat ± 0.17 sys in minimum bias and central (0-20%) d+Au 
collisions, respectively. 
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Cross sections at large transverse momentum (px) for 
products of "hard" point-like processes (e.g., inclusive 
hadrons, jets, direct 7's, and heavy flavor) in high en- 
ergy hadron collisions are well described in perturba- 
tive Quantum Chromodynamics (pQCD) As a re- 
sult they are considered to be well-calibrated probes of 
small-distance QCD phenomena. When going from p + p 
to p(or d)+A and A+A collisions, deviations from cross 
section scaling with respect to the number of binary N+N 
collisions provide information on cold nuclear matter ef- 
fects such as initial state energy loss 0, shadowing Q 
and hot nuclear matter effects such as in-medium energy 
loss [4|, increasing importance of production via recom- 
bination [f|, and modifications to the QCD vacuum [f|. 

Production of to mesons at high-px is especially in- 
teresting. The to and tt° are vector and pseudoscalar 
mesons, respectively. The to/n ratio carries information 
about probabilities of corresponding spin states to be pro- 
duced in hadronization. Furthermore, to mesons can be 
used as a probe of the nuclear medium since a significant 
fraction of them produced in A+A or p(d)+A collisions 
will decay inside the produced medium (ct = 23.8 fm) 
possibly leading to changes of the mass and/or width of 
the to with respect to their values in vacuum. 

The PHENIX Q experiment at RHIC has a unique 
capability to measure both neutral and charged prod- 
ucts of A+A and p + p collisions at very high event 
rates. The two central arm spectrometers each cover 90° 
in azimuth and ±0.35 in pseudorapidity. The Electro- 
magnetic Calorimeter (EMCal), with energy resolution 
a/E = 8.1%/y/E(GeV) © 2.1%, is used to reconstruct 
7's and 7r°'s. For charged particle reconstruction two 
layers of Pad Chambers provide 3D pattern recognition 
and fake track rejection, and a Drift Chamber gives mo- 
mentum resolution a/p T = 0.7%y/pr(GeV/c) © 1.1%. 
Beam-Beam counters (BBC) were used to provide mini- 
mum bias trigger and to determine the collision vertex. 
The minimum bias trigger cross sections measured by the 
BBC are 23.0 + 2.2 mb in p + p 8] and 1.99 + 0.10 b in 
d+Au [9) collisions. In d+Au collisions the BBC were 
also used to separate events into centrality classes as ex- 
plained in [10| | . High px online triggers are implemented 
by adding together amplitudes in 4 x 4 adjacent EMCal 
towers and comparing them to a threshold of 1.4 GeV in 
p + p and 2.4 GeV in d+Au. The trigger could be fired 
by one or more photons coming from to decay final states 
including 37 or 27r27. 

The data was collected by the PHENIX experiment 
during RHIC Run3. After selecting good runs and cut- 
ting on the collision vertex (|z| < 30cm) we analyzed 
approximately 4.6 x 10 7 and 2.1 x 10 7 high px trigger 
events, corresponding to a total integrated luminosity of 
0.22 pb _1 and 1.5 nb -1 in the p + p and d+Au collision 
systems, respectively. 

Although w-mesons are relatively abundant in high en- 
ergy hadronic collisions (to/jr w 1 at high px), their 



measurement is challenging due to the multi-particle final 
states of the main decay channels and the combinatorial 
background associated with their reconstruction. 

The procedure used to measure w — > 7r°7r + 7r _ is the 
same as used to measure rj — > 7r°7r + 7r~ 10]. In case 
of to one shall account for a wider peak and differ- 
ent phase space density of the 3-body decay 
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The analysis procedure for the photonic decay mode, 
to — ► 7r°7 — * 37, is very similar to earlier PHENIX mea- 
surements of other mesons with photonic decay modes, 
7T . 77 — > 27 [j.[H.[l3|. For both to decay modes studied 
the first step is to reconstruct ir° candidates by combining 
photon pairs and applying a px dependent cut around the 
mass of the ir°. The R.M.S. of the ir° peak varies with px 
from 8MeV to 13MeV. Candidates (which include com- 
binatorial background) are combined with a third photon 
for to — > 7r°7 or with two unidentified charged tracks (as- 
sumed to be 7r-mesons) for to — > 7r 7r + 7r~. Raw yields are 
extracted by fitting the px slices of the invariant mass 
distribution as in the insert panel in Fig. [1] The signal 
to background ratio (S/B) of 7r°7r + 7r~ decay channel is 
1:7 at low px and grows to 1:1 with px in p + p collisions. 
In d+Au collisions it starts at 1:20 and grows with px 
to 1:2. For the 7r°7 channel S/B is 1.5-2 times worse. 
Corrections for acceptance, trigger efficiency, and anal- 
ysis cuts are described in detail in [l(|. Further details 
about the analysis procedures can be found in [l5j |. 

We classify systematic error sources as type A (point- 
to-point uncorrelated, which can move each point inde- 
pendently), type B (point-to-point correlated, which can 
move points coherently, but potentially by different rela- 
tive amounts), and type C (global, which move all points 
by the same relative amount). These errors are summa- 
rized in Table [J for the different decay modes and colli- 
sion systems. Major contributors include signal extrac- 
tion (type A), online high-px trigger efficiency (type B) 
and the total cross-section measurement (type C). The 
uncertainty on the signal extraction, which is the dom- 
inant source of systematic error, is relatively large due 
to the fact that correlations in the triggered sample ren- 
dered background subtraction via event mixing impossi- 
ble, and therefore the unknown shape of the underlying 
background had to be accounted for in the peak fit. This 
error is estimated based on variation of analysis cuts and 
fitting procedures (loj . 

Figure [T] shows the invariant cross sections for to pro- 
duction in p + p, minimum bias d+Au, and central (0- 
20%) d+Au collisions at ^/snn = 200 GeV, as a function 
of px- The px range (2.5 < px < 9.25GeV/c) is limited 
by statistics at high px, and by decreasing detector ac- 
ceptance and trigger efficiency at low px- The results for 
the two decay modes, which involve different kinemat- 
ics, detector acceptance, and efficiency corrections, agree 
very well. 

Figure shows the results for the ratio of vector to 
pseudoscalar meson production (00/7T ) in d+Au and p + 
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TABLE I: Systematic errors for the u production cross sec- 
tion for the two decay channels and the two collision systems 
analyzed. Error types are described in the text. Values with 
a range indicate variation of the systematic error over the px 
range of the measurement. 
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FIG. 1: Invariant cross section of uj production in p + p 
and d+Au collisions at ^Jsnn = 200 GeV measured in cu — > 
7r°7r + 7r~ and u — > 7r°7 decay channels. Bars and boxes rep- 
resent statistical and systematic errors, respectively. Fits to 
invariant mass distributions of (top) and 7T°7r + 7r (bot- 
tom) are shown in the insert. 



p collisions at y?SNN = 200 GeV. For the denominator 
we use inclusive ir° yields measured by PHENIX [H, [l(| . 
Ratios in both systems are consistent with unity over all 
measured pp. Fits to a constant yield oj/ir° = 0.85 ± 
05 stat ± o ogsys an d .94 ± 0.08 stat ± 0.12 s y s in p + p 
and d+Au collisions, respectively. Fits assuming linear 
Pt dependence have slopes consistent with zero. At high 
Pt the PYTHIA [13] prediction for this ratio in p + p 
collisions at y/s = 200 GeV is consistent with, but slightly 
higher than, the measurement. 
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FIG. 2: Measured u>/ir vs pt in (upper panel) d+Au and 
(lower panel) p + p collisions at ^/snn = 200 GeV. Straight 
lines show fits to a constant for each collision system. The 
boxes at the left edge of the constant fit lines show the sys- 
tematic error on the data averaged over px- Dashed lines 
show values within la of the best linear fits to the data. The 
PYTHIA prediction [13] for p + p collisions at Vs = 200 GeV 
is shown with a solid curve in the bottom panel. 



The R-806 experiment at ISR measured u>/ir° in p + p 
collisions at y/s — 62 GeV [ll| and found this ratio 
to be 0.87 ± 0.17 over 3.5 < p T < 7GeV/c. The 
E706 experiment measured ui/tt° in 7r~+Be collisions at 
yfs^N = 31 GeV (l9j and found values consistent with the 
results presented in this paper. The uj/it ratios mea- 
sured in hadronic interactions by three experiments at 
three different energies between 31 and 200 GeV, are the 
same within the errors. 

Several LEP experiments 2(], 21, 22 1 have measured u 
production in e + + e~ collisions at \fs = 91.2 GeV and 
presented results as a function of x p = 2p UJ /y/s. This 
is not a well-defined quantity in hadronic interactions. 
However, one can compare the tu/ir ratio at large values 
of pt and x p . Following the procedure for handling LEP 
data described in [l(| one finds that for x p > 0.5, the 
largest value for which statistically significant LEP data 
is available, the ratio has grown to approximately 0.7, 
close to the measurements in hadronic interactions. 

Figure [3] shows the nuclear modification factor i?^ A , 
denned as the ratio of the u> meson yields in d+Au in- 
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teractions and p + p interactions scaled by the num- 
ber of binary collisions in d+Au, for minimum bias and 
central (0-20%) d+Au events. Precise definition of the 
RdA and procedure to determine centrality in <j+Au is 
given in sections IV.B and III.C in [k]]. We find R% A 
to be 1.03 ± 0.12 stat ± 0.21 sys for minimum bias and 
0.83 ± 0.21 stat ± 0.17 sys for central events, independent 
of pt . The RdA for two other neutral mesons (ir° and 77) 
measured by PHENIX 0, El are also shown in Fig. [3] 
In all cases PHENIX observes that RdA is close to one 
for pt > 2GeV/c and flat out to the highest pr- Simi- 
lar behavior is seen in preliminary analysis of K and <p 
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FIG. 4: Reconstructed u> mass measured in the 7r°7r + 7r _ decay 
channel vs pt in (upper panel) d+Au and (lower panel) p + p 
collisions at ^/snn = 200 GeV. Error bars show statistical 
errors. Straight lines show fits of the data to a constant. 
Dashed lines show values within la of the best linear fit to 
the data. Systematic error on the fit value is shown with the 
box. PDG values for uj meson mass (m) and width (r) are 
shown with arrows on the left I27I1. 
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FIG. 3: Measured RdA vs pt for neutral mesons in d+Au col- 
lisions at y / sjvjv = 200 GeV for (upper panel) 0%-20% central 
and (lower panel) minimum bias. Values for 7r°'s and ?7's are 
from [id 1 16j| - For reference a line is plotted at RdA=l- The 
scaling systematic error is shown as a box on the left. 



Recent publications suggest that modifications to the u 
mass can be observed even in cold matter by studying not 
only the electron decay channel [24],[25j], but also hadronic 
channels [26j]. For the hadronic decay modes presented 
here PHENIX lacks acceptance at low pt where the effect 
is expected to be the most prominent. However, we do 
have excellent mass resolution (20-25 MeV) for the mixed 
neutral-charged particle decay mode. Figure [4] shows ex- 
tracted values for the u> mass as a function of pr- I n the 
Pt range of the measurement we observe no modification 



In summary we have presented the first measure- 
ment of to production in p + p and d+Au collisions at 
^snn — 200 GeV. The production cross section is mea- 
sured in two different decay modes with consistent re- 
sults. The ui/tt ratio in p + p collisions is found to be 
0.85 ± 0.05 stat ± 0.09 sys and 0.94 ± 0.08 stat ± 0.12 s y s in 
d+Au. over the measured range. This agrees with 
previous measurements in hadronic collisions at lower 
y/s. The nuclear modification factor for u> production in 
d+Au collisions is consistent with 1 and pt independent 
for pt > 2GeV/c consistent with other meson measure- 
ments. No modifications to the u> mass were observed in 
p + p or d+Au collisions. The u) meson is also interesting 
as a probe of the hot nuclear medium created in A+A 
collisions. This measurement will serve as an excellent 
baseline for measurements of u> production in A+A in 
various decay channels. 
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